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Abstract 

' For scalar reaction-diffusion in one space dimension, it is known for a long time that 

fronts move with an exponentially small speed for potentials with several distinct mini- 
mizers. The purpose of this paper is to provide a similar result in the case of systems. 
Our method relies on a careful study of the evolution of localized energy. This approach 
Ph ' has the advantage to relax the preparedness assumptions on the initial datum. 
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1 Introduction 

1.1 Potentials with wells of equal depth 



The purpose of this paper is to investigate the behavior of solutions v of the reaction-diffusion 
equation of gradient type 

y^-y^^ = -^V{v). (1) 

in 

, The function v denotes here a function of the space variable x E M and the time variable 

<^ I t > and takes values in some euclidean space M^, so that ([1]) is a system of k scalar partial 

■ differential equations. Equation ([T]) actually corresponds to the gradient-flow of the energy 

functional £ which is defined for a function n : M i-^ M'^ by the formula 

X: £^u)= [ e{u)= [ \^ + V{u). (2) 

The function V, usually termed the potential, is assumed to be a smooth function from M.^ 
to M, tending to infinity at infinity, so that it is bounded from below. 

Simple solutions to equation ([1]) are provided by the stationary ones, that is time-independent 
solutions of the form v{x, t) = u{x), where the profile n : R i-^ M'^ is a solution of the ordinary 
differential equation 

-u^x = -yV{u). (3) 

Among those solutions, the simplest ones are the constant functions v{x,t) = a, where a 
is a critical point of the potential V, for instance a minimizer. Another interesting type of 
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solutions to ([3]) is provided by those tending, as x — > ±00, to critical points of the potential 
V: in this case conservation of energy for ([3]) implies that V{u{+oo)) = V{u( — oo)). The 
central assumption on the potential V in this paper is that it possesses a finite number of, 
and at least two, distinct minimizers. A canonical example in the scalar case A; = 1 is given 
by the function 

V(u) = ^. (4) 

whose minimizers are +1 and —1. 

This paper is devoted to the analysis of the evolution in time of initial data which connect 
two distinct minimizers of the potential. Such maps u, from M to M'^, whose limits at ±00 
are distinct minimizers of V are usually termed fronts. If they are moreover solution to the 
ordinary differential equation we will call them stationary fronts, so that a stationary 
front is a heretoclinico solution to ([3]). For instance, in example ([4]), stationary fronts are 
necessarily of the form 

(X — c\ 
~~\/2~ ) ' ^^'^ 

for some c S M. In that case, (resp. w^), which is often referred to as the kink (resp. 
anti-kink) solution, connects —1 to +1 (resp. +1 to —1). The dynamics of fronts and their 
eventual convergence to stationary fronts, which are attractors of the dynamics, is actually 
a central topic in the study of reaction-diffusion equations of gradient-type: in our short 
historical section below, we will review some of the works related to this question. 

To be more specific, our assumptions on the potential V can be formulated as follows. We 
assume that V is smooth and satisfies the conditions: 

(Hi) inf y = and the set of minimizers S = {y G R'', V{y) = 0} 

is a finite set, with at least two distinct elements, that is 

S = {ai, aj, q>2, cTi G M^ = 1, q. (6) 



(H2) The matrix V^F(ai) is positive definite at each point (Ji of S, in other words, if 
denotes its smallest eigenvalue, then A~ > 0. We denote by its largest eigenvalue. 

(II3) There exists constants oq > and Rq > such that 

yVV{y) > ao|?/|^ if \y\ > Ro- 

A potential V which fulfills conditions (Hi), (H2) and (H3) will be termed throughout a 
non-degenerate multiple-well potential with equal depths. A canonical example is given by 
gl, for which S = {+1, -1}. 

The main assumption in this paper on the initial datum v^{-) = v{-, 0) is that its energy is 
finite. More precisely, given an arbitrary constant Mq > 0, we assume that 

(Ho) £{v°) < Mo < +00. 

^Actually homoclinic solutions, whenever they exist, could be considered as well. 
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In particular, in view of the classical energy identity 



£{v{;T2))+ r 



dv 



we have, Vt > 0, 



{x,t)dxdt = £{v{-,Ti)) V0<Ti<r2 



£ivi-,t)) < Mo. 



(7) 



(8) 



This implies in particular that for every given t > 0, we have V{v{x,t)) — > as — > oo. 
It is then quite straightforward to deduce from assumption (Hq), (Hi), (II2) as well as the 
energy identity ([7]), that v{x,t) (X± as x ^ ±00 , where a-t S S does not depend on t. 
In other words our assumptions imply that the map v{-,t) is a front for all times t > if 



1.2 Front sets 

One of our aims is to localize the evolution in time of the region where the function v{-,t) 
jumps from one minimizer of 1^ to a second one. This will allow us to follow the evolution of 
the front. To that purpose, we fix |j.o > sufficiently small so that, for i = 1, . . . , we have 

for all i ^ j in {1, • • • ,q} and 

^Xrid < vV(y) < 2A+Id (9) 

for all i G {1, • • • , q} and y G i?(o',j, \Xq). We then define, for a map n : M 1— > M'^, the set 

V{u) = {x G M, dist(n(x), S) > ^o}• (10) 
In the context of equation ([1]) we set moreover 

Vit)=V{v{.,t)). (11) 

The evolution of the set "D^t) is the main focus of our paper. 

For a given map u, T){u) is related to the set where the energy of u concentrates, in view 
of the following: 

Lemma 1. There exists a constant r/o > 0, depending only |j.o > and V, such that, if I is 
an interval ofM. of length \I\ > 1, and u is a -valued function on I satisfying 

eiu) < (12) 



then 

dist(u(3;), S) < iJ-o for all x e I, (13) 

or equivalently 

V{u) nl = $. (14) 
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This kind of result is usually called a clearing-out Lemma in the literature. It shows that 
if the energy is sufficiently small in some place, then there are no front located there, or 
equivalently that where fronts are present, energy needs to concentrat^. Therefore, fronts 
are among energy concentration intervals, and energy is a good object to track fronts. We 
will explicitly assume in the sequel Mq > r/o- 

An immediate consequence of Lemma [T] yields: 

Corollary 1. Assume that the map u satisfies £{u) < Mq. There exists £ points xi, ...,Xi in 
'D{u) such that 

C Uf=i[xi- + (15) 
with a bound I < — on the number of points. 



1.3 Slow motion of concentration sets 

The first and main result of this paper is as follows. Assume Mq > rjQ, and set 

Mo 

ao = 32 . 

Vo 

Theorem 1. Assume that the potential V satisfies assumptions {Hi), {H2) and (-H3), that the 
initial datum vq satisfies the energy bound (Hq). There exists a constant Kq > depending 
only on the potential V and on Mq such that if R> ctQ, 

V{t) cV{0) + [-R,R] (16) 

provided 

rV f r 



Theorem [J expresses the fact the motion of the front set is slow when considered at suffi- 
ciently large scale R. Indeed, its maximal average speed should not exceed 



c{R)=KlR-^e^vl^-^^ 



For large R this speed is exponentially small. As a matter of fact, our proofs rely on an 
asymptotic expansion for large values of R. 

Since R represents the typical length scale which is considered, if one wishes to work on 
fixed domains, it is sometimes interesting to introduce a small parameter e > 0, and consider 
the more general form of the equation ([1]) given by 

dtVe- dxxVe = --^yViVe). (17) 

Notice that if t; is a solution to ([1]), then the map given by 

Ve{x,t) = v{ex,e'^t) (18) 



The converse is of course not true in general for arbitrary maps, think of smaU osciUations. 
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is a solution to ()17l) . In this setting the statement of Theorem [T] may easily be translated if 
we replace equation ([1]) by equation (fT7|) and the map v by the map Ve- The energy density 
es{u) and the energy functional <fe(n) are given respectively by 

e,(u) = Jn^ + ly(^) ^nd f,(n) = / e,(n). (19) 

The energy bound (Hq) is translated into 

(Ho") Seiv'^) < Mo < +00. 

We may define accordingly the front set T>{t) = T>{vs{-,t)), so that if the initial datum satisfies 
the initial bound (Hq)", this set is of size of order e, in view of Corollary [T] and (fTS]) . and 
shrinks, as e converges to zero, to a finite set, which is sometimes termed the defect set. 
Considering R > fixed but letting e > vary, the statement of Theorem [T] turns then out 
to be equivalent to the following statement. 

Theorem Ibis. Assume that the potential V satisfies assumptions (Hi), {H2) and {H^), let 
e > be given and consider a solution to (jl7p . Assume that the initial datum v^{-) = 
Ue(-,0) satisfies the energy hound (Hq). There exists a constant Kq > depending only on 
the potential V and on Mq such that if R> oiqE, then 

V{t) (ZV{<d) + [-R,R] (20) 

provided 

rV f r 



< t < — exp 



KqJ \Ko£ 

As we will recall in our short historical survey on the topic below, slow motion of fronts 
has a long history in the mathematical literature, and it is mainly described considering the 
form ([1]) and assuming e is asymptotically small. 

Combining Theorem Ibis with dissipation estimates off the front set (see Section U] below), 
we obtain 

Theorem 2. Assume that V and are as in Theorem Ibis. There exist a constant Ki > 
depending only on V such that if xq gE. and R > oqE satisfy 

[xo - 2R, xo + 2R] n P(0) = 0, 

then 



exp —rr-o +7^ exp 



R 



Kisy ' R^ Kis 

pointwise on [xq — ^R,xo + ^R] x [e^, l"^)^ exp(-^^)]. 



(21) 



Using a Gronwall type argument, we may then prove that the fiow drives the solution close 
to a chain of stationary fronts. More precisely, we have 
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Theorem 3. Assume that V and are as in Theorem Ibis. Given R > ao£ there exists a 
relaxation time 

rV f r 



< T < ( — ) exp 

foil 

and Mq, a length scale of order R 



at which Vi;{-,T) possesses the following structure: There exist K2 > depending only on V 
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M, 



_,_o R R 
2 'w — < r < — . 
K2- - K2 



a collection of points {ojjjgj mM, and a corresponding collection of functions defined 
on [— ^, ^] with values into M'^ such that 

1. < aj < 

2. aj G V{T) V i G J, 

3. dist(aj,P(0)) < R, V j G J, 
4- dist(aj, Oj) > 4r y i ^ j £ J, 

5. Each Uj is a solution to the stationary equation ([3]) with zero discrepancy: 



|2 



-d.xU, = -VV{Uj), i{Uj) = - V{Uj) = 



6. We have the estimate 

\\v,{.,T)-U,(^-^) \\+e\\d^(^Ve{;T)-U,{^-^)) II <i^2exp (-;^) > 

in L°° ( [oj — r, aj + r] ) , for each j G J, 

7. If I is an interval disjoint from Ujg j[aj — r, aj + r], we have 

\\ve{.,T)-a,\\+e\\d^Ve{;T)\\ < K2exp 

in L°°{I), for some cij G S. 

Notice that the functions Uj are defined on the interval [— ^, which grows as e tends to 
if R is kept fixed, to cover the whole of M. In particular if one considers a family of solutions 
(''^e)o<£<i to pT|) satisfying the energy bound (Hg) and the corresponding family Uj = Uj 
obtained thanks to Theorem [3l then a straightforward compactness argument shows that, up 
to a subsequence e„ ^ 0, we have the convergence 

^ ' as n ^ 00, (22) 

in C^{K) for any A; G M and any compact interval K of M, where the limiting map Uj is 
defined on the whole of M and is a stationary front, that is a finite energy solution to ([3]). 
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Loosely speaking^, one may rephrase Theorem [3] stating that, after some suitable time, the 
solution enters an O (exp(— -^)) neighborhood of glued together stationary fronts. 

The next step in the analysis would be to derive a precise motion law for the fronts, or more 
precisely, in view of Theorem [3l of the points aj . It is in particular of interest to determine 
whether the interactions between them is repulsive or attractive. Notice however that the 
points aj are only defined so far up to an 0{R) term. This introduces an additional difficulty 
which is easily removed in the scalar case imposing some preferred value for C/j(0), in the 
vectorial case the situation is more subtle. We are inclined to believe that such results would 
involve more restrictive assumptions on the potential V than the ones which we have used 
so far, which are rather mild and involve only its behavior near its zeroes. In particular, 
more should be known or required regarding the zoology of finite energy stationary front 
(concerning for instance their spectral properties). 

In the scalar case, it is a standard exercise to integrate equation ([3]) and to determine the 
set of finite energy solutions (for instance, the solution is given by formula ([5]) in the case 
the potential is (j3])). Furthermore in that case, fixing the discrepancy to zero as it is done in 
statement 5 of Theorem [3] insures that all the solutions Uj are true stationary front^ (that is 
they coincide with finite energy globally defined solutions of dS])). In a forthcoming paper, we 
will show how these information combined with the local energy identity which is the central 
tool of the present paper allows to recover the motion law in the Allen-Cahn case (see the 
historical notes below), where the fronts attract, and also can be extended to three wells or 
more, where the interactions may be attractive or repulsive, depending on the nature of the 
fronts. 

An ultimate goal would be of course to obtain similar results in the case of systems, when 
appropriate assumptions are made on the potential. In that case, we believe that the in- 
teraction between two fronts is governed by the behavior as x — > ±00 of the corresponding 
finite energy stationary fronts. In order to state a conjecture which specifies the magnitude 
of the interaction, consider first a given finite energy stationary solution [/ to ([3]), and set 
(T^ = lim U{x). In view of the form of the potential near its zeroes, it can be shown that the 

X— >itoo 

solution converges to cr"^ (resp. <s^) with an asymptotic direction, that is that the following 
limits exist 

CU+ = lim ^^-^ ^ fresp. tu" = lim ^^— \ , 

x^+oo|[/(x) — a+l \ x^-<x\U{x) — J 

so that I a)"*" I = |cu~| = 1. Next, we go back to Theorem [3] and consider two consecutive 
"fronts", say Ui and U2 given by its statement. We have in particular a~^{Ui) = g~{U2) = c. 
Let A > denote the smallest eigenvalue of the matrix V^y(cT), let F denote the correspond- 
ing eigenspace, and 11 the projection onto F. Our conjecture is that the interaction between 
the two fronts is given by 

(C(ntuJ^, IltUg + 0(1)) exp ^— — , as e ^ 0, 
^Loosely in particular because the f/j's are not globally defined ! 

*In particular, for the potential given by Q , Theorem |3] allows to recover the convergence result of Fife 
and McLeod [11] without any reference to the maximum principle (this was already achieved in Gallay and 
Risler [13]), and with a quantitative error estimate (we do not rely on a compactness argument). 
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where the constant C > depends on the profiles Ui and U2 and d denotes the distance 
between the two frontsU 



A short historical review. General reaction-diffusion of gradient type ()17p have been 
widely introduced and used as models in numerous branches of science, for instance in physics, 
chemistry (in particular combustion theory), or biology, among many others. The theory of 
fronts has received extensive mathematical study, in particular the theory is highly developed 
in the scalar case. For instance, existence and uniqueness (up to translation) of traveling 
front^ has been established in the scalar case. It has also been shown that for arbitrary 
initial data connecting the local minimizers, as time tends to infinity solutions converge 
towards such fronts (see e.g. the seminal work of Fife and McLeod jllj). More recently, 
part of the analysis of convergence towards traveling fronts as time tends to infinity has been 
extended to the case of systems by Risler (2Tl[22l[23] and Gallay and Risler[13]. The methods 
used in these works rely on energy estimates and compactness arguments. 

When a front connects two local minimizers with the same potential energy, it is stationary, 
and therefore is a solution of (l3|). As proved by Risler, as time tends to infinity solutions 
for "arbitrary" initial data eventually converge to slowly repulsing chains of such stationary 
fronts. The next step of the analysis, when the initial data possesses several fronts, is to follow 
carefully the evolution of the various fronts from the initial time, and not only asymptotically, 
and possibly to estimate their speed. On a heuristic level, the speed of the fronts can be seen 
as the effect of the (small) interaction between them: This interaction might be attractive or 
repulsive depending on the nature of the potential. 

For scalar two-wells potentials, which are often referred to as Allen-Cahn potentials, this 
program was first completed in the celebrated works of Carr and Pego, and Fusco and Hale 
[SI [71 [12], which provided the first rigorous mathematical derivation of very slow motion, 
and even derived a precise motion law for the evolution of fronts. In their result, the initial 
data is very constrained, since it is supposed to be close to optimally glued together fronts. 
Their method relies on a careful analysis of the motion near these special solutions through 
a thorough study of the linearized operator near the stationary solution. The fact that the 
kernel for such solutions on the line contains only the space derivative of the solutions is 
crucial there (the proof of this latest statement provided in [6| relies heavily on the fact that 
the solution is scalar). This type of spectral methods were later applied successfully on related 
problems, for instance the Cahn-Hilliard equation (see e.g. [Il[2])- Other interesting papers 
based on that kind of ideas (sometimes termed the invariant manifold method or geometric 
method) are [9l [HI [T2| [20l [TO] . In particular in [1^ , Ei was able to handle the interaction of 
a kink and an anti-kink in the vectorial case. Keller, Rubinstein and Sternberg [16] made 
related contributions in a similar direction. 

At least two other methods have been applied successfully in the scalar case. Firstly, the 
method of sub-solutions and super-solutions turns out to be extremely powerful and allowed to 
handle larger classes of initial data (see e.g. |1H[8]). There is little hope however to extend this 
method to systems, since comparison principles do not hold in general for systems. Another 
direction is given by the global energy approach due to Bronsard and Kohn [5]. Using the 
energy identity ([7|), they were able to prove that for initial data sufficiently close to glued 

^As explained above, defining the distance between two fronts supposes first to localise them by fixing an 
anchor point for each of them, the constant C depends of course on the definition of those anchors. 
®that is, solutions with a constant profile which are translated with constant speed 
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front solutions, the fronts have a speed slower than 0{e^), for every /c G N. The closeness to 
the glued solutions is expressed in their paper in terms of energy estimates and assumptions 
which are reminiscent of concepts of Gamma-convergence. Grant [Hj improved the method 
to obtain an exponentially small upper bound of the form O (exp(— |)) imposing however 
stronger conditions on the initial data. The method was extended to functionals with higher 
order derivatives in [2l]. Finally, in [19] through a more abstract setting of the problem. Otto 
and Reznikoff were able to recover some of the results in [8] through global energy methods. 

At this stage it is worthwhile to emphasize that all mentioned results use more or less the 
properties of the solutions to ^ (in the case of Allen-Cahn, these are unique up to transla- 
tions and sign change). 

The aim of our present paper is to extend the analysis to the case of systems, to relax the 
assumptions on the preparedness of the initial data, and also to handle possible annihilation 
of fronts. Our approach bears many analogies with the energy method of Bronsard, Kohn 
and Grant, however, instead of using the global energy identity we use a local version of it, 
which is combined with parabolic estimates away from the fronts. Our ideas are partially 
borrowed from our earlier work on the motion of vortices in the two-dimensional parabolic 
Ginzburg-Landau equation [31 [4] as well as from earlier works on the topic by Lin, Jerrard 
and Soner [17^ I15j . A new technical difficulty which occurs in the present paper is that we 
are able to handle, for the evolution of localized energy density, test functions which are 
affine (and therefore change sign) near the front, whereas (positive) quadratic test functions 
were extremely useful in the context of the Ginzburg-Landau vortices. In contrast with the 
results obtained so far in the scalar case, our results do not rely on the properties nor the 
existence of solutions to ([3]) (as a matter of fact, the number and the properties of solutions 
to ([3]) might be much more involved in the multi-dimensional case than in the scalar case). 
Roughly speaking the heart of our method is precisely to avoid regions where the solution 
becomes close to solution to ([3]). 



1.4 Elements in the proofs 

For the proofs of the main results, we will work with the parameter e (equation ([T]) being a 
special case for the value e = 1). The most difficult part corresponds to the case where e is 
small with respect to R. 

In order to analyze the evolution in time of the concentration sets of the energy, we invoke 
the localized version of ([7]), which writes, for a smooth test function x with compact support 
in M 



dt 



I ea.v).(.)<J. = -/ .a,„.V,,.Vx<i.. (23) 

JRy.{t} JRx{t} JRx{t} 

A few integration by parts yield the classical formula (|24p below. 

Lemma 2. Let x be a smooth function with compact support on R. Then, we have the 
identity 

— I X{x) ee{ve)dx = - i ex{x)\dtVe\'^dx + J^s{t,X,Ve), (24) 
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where, the term Ts, is given by 



x{t} 



■ V(vs 



e 



X dx. 



(25) 



The first term on the right hand side of identity (j24p stands for local dissipation, whereas 
the second might be interpreted as a flux. The quantity 



ViVs 



(26) 



is sometimes referred to as the discrepancy term in the literature. For solutions of the 
ordinary differential equation 



+ \vV{ue) = on /, 



(27) 



for some interval /, the discrepancy is constant, and in particular it vanishes if the interval 
is the whole of M and the solution is a heteroclinic between two of the wells. A general idea, 
which is underlying our analysis, is the fact that front sets in the parabolic equations relax 
quickly to solutions to ()27p . so that we may take advantage of some properties of the ode 
()27p . in particular in order to estimate the driving force J-'s- 

We will make use of formula (I24p for a rather specific choice of test functions X: namely test 
functions which are affine near the defect set, so that x vanishes there, and one has merely, 
in view of Lemma O to estimate the discrepancy term off the defect set. An important step 
in our proofs is to show that if the solution is locally close to one of the minimizers cTj at a 
given time, it remains so for a short while, and the equation has strong smoothing properties 
on the corresponding space-time region. 

We have 

Proposition 1. Assume that Vg is a solution to (jl7p verifying assumption (Hq). Let t >0, 
a; G M and r > oqE be such that 



Then, 



'D{t) r\[x — r,x + r] 



V{s) n [x - r/2,x + r/2] 



(28) 
(29) 



for every t < s <t + ^r^. Moreover, for such s. 



x+r/2 \ + 

ee{ve{y,s))dy<Mo{l + A^) 

x-r/2 



exp(-%(s-t)) + 2 
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i - t , / A,- 1 r , 



where \- and Xf are defined in {H2) for the corresponding potential well. 



(30) 



Finally, another important tool in our proofs in an elementary covering argument, which 
we introduced in an earlier work on the motion of Ginzburg-Landau vortices and that we will 
recall in subsection [5j Roughly speaking, given a small number k, this covering results states 
that given I points in a metric space and a typical length 5, one may find a length scale 5 
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of the same order as 5, such that all the points are included in balls of radius k6, and such 
that mutual distances between the balls is at least k~^6 (the points are said to be K-confined 
inside these balls). 

The proof of Theorem Ibis is then a combination of applications of Lemma [21 Proposition 
[U the clearing-out Lemma [H and the K-confinement Lemma. It is based on an iteration 
argument, which will be completed in a finite, bounded by above, number of steps. The main 
parameter which is involved in the iteration process is 6, which appears in the n- confinement 
technique, and which we let it shrink. More precisely, we set Sn+i = SnK^^', and each step on 
the iteration is defined thank to a suitable stopping time. 



2 Properties of the front set 

We describe in this section some classical and mostly elementary properties of the front set. 
More precisely, we present the proofs of Lemma [T] and Corollary [TJ 

Proof of Lemma [H Let a G / be such that [a, a + 1] C /. Since by assumption \I\ > 1, 
the set of such points is non-empty. In view of the definition of the energy, V{u) < e{u), and 
therefore if the map u satisfies (I12p then 

V{u{x))dx < / V{u{x))dx < rjQ. 

[a,a+l] Jl 

We deduce from the mean- value theorem that there exists some point sq G [a, a -|- 1] such 
that 

V{u{so)) < m- (2.1) 

We first choose the constant r/o small enough so that V{y) < rjQ for y G implies that 
y G B{ai, \Xo/2) for some z G {1, • • • ,q}. In view of (j2.ip . there exists therefore i £ {!,■ ■ ■ ,q} 
such that 

Hso) - ai\ < ^. (2.2) 
On the other hand, by integration we have for any s G [a, a + 1] 

\uis)-uiso)\< [ \u\<{[ |n|2)^ < (277o)i (2.3) 

J[s,so] -'[s,so] 

We impose additionally that (2ryo)2 < so that combining (|2.2p and (|2.3p we obtain, for 
any s £ [a,a + 1], 

\u{s) - Gil < ^. (2.4) 

Since this relation holds for any a € I such that [a, a + 1] C I, inequality (|2.4p actually holds 
for any s G /, and the proof is complete. □ 

Proof of Corollary [H We consider the covering of R given by R = U /„, where In denotes 
the interval [n, n -|- 1]. Consider the subset T of Z defined by 

I = {n£Z, V{u)nln + 0}. 
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In view of Lemma [T] we have fj e{u) > r/o for any n £ I, so that we deduce 

Mo> f e{u) >Y. I ^(^) ^ (W^)^O' 



neX' 



so that < For n £Z choose Xn G ^(^i) H /„, so that In C [xn — 1, x,i + 1]. We have 

V{u) C J [x„ - 1, x„ + 1] . 

»i£X nSX 

□ 

3 A first upper bound for the velocity of the front set 

In this section we provide a first estimate concerning the speed of the front set. It relies on 
the clearing-out lemma as well as on some maximum principle for the quantity | f g — CTj p . The 
main result of this section is the following 

Proposition 3.1. Assume that Ve is a solution to (I17p verifying assumption {Hq). Lett > 0, 
X and r > a^e be such that 

V{t)n[x -r,x + r] =(/}. (3.1) 

Then, 

3 3 

V{s)n[x--r,x + -r] = (3.2) 

for every t < s < t + ^r^. 

The starting point is a straightforward consequence of formula (j23p . often termed semi- 
decreasing property. 

Lemma 3.1. Let x be a smooth function with compact support in M, then 

\f e\dtVeW^ + ^ [ ee{v,)x^ < A\\x\\l^ [ e,{ve). (3.3) 

^ jRxH} JRxH} -'suppy 



2 

In particular, 



di 



I ee{ve)x^{x) dx<A llxllioo £M) < 4Mo||x||ioc, (3.4) 
JRxH} 



<{t} 

so that forO<t<t + At 

[ eeiveix,t + At))x\x)dx< [ eeiveix,t))xHx)dx + 4Mo\\x\\looAt, (3.5) 

JR JR 
whereas in the opposite direction 



ee{ve{x,t + At))x {x)dx> I ee{ve{x,t))x {x)dx 

■t+At 



(/ e\dtVe\\^)-4M4x\\looAt. 

Jt JR 



(3.6) 
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This property yields 
Lemma 3.2. Under the assumptions o f Proposition [0|, there exist 

Xj E [x — r,x — ^r] and ^ ~^ x + r] 

such that 

{xi,x/}nP(s) = 

for every t < s < t + Oq ^r^. 

In other words, the two time segments {xj} x [t,t + a^^r'^] and {xf} x [t,t + ag^r^] are 
disjoint from the space-time front set of v^- 

Proof. By the pigeonhole principle there exists Xj E [x — r, x — |r] such that 

° ee{ve{y,t))dy < —. 

2r 2 

atQ 

Indeed, the total energy of Ve is bounded by Mq and an interval of length r/4 is ao/16 times 
larger than an interval of length 4r /ao • 

Let X be a smooth non-negative cut-off function such that x = 1 on [xj — e, Xj + e], x = 
outside [x — 2r/ao,Xj + 2r/ao], and |x'| < aor~^. Applying (13. 5p we obtain 

ee{ve{y, s)) dy<'^ + 4Moalr-\s - t) < % 



provided 



One obtains similarly Xj and the conclusion follows from Lemma [TJ □ 

By continuity, it follows from previous lemma that under the assumptions of Proposition 
13.11 there exists (7j E S such that takes values in -B((Ti, |Xo) on the parabolic boundary 

dpk = [xi,Xf] X {t} \J{xi} x[t,t + ao V^] |J{x/} x[t,t + Qq V^] 

of the cylinder 

A = [xi,x/] X [t,t + OQ^r'^]. 

We are now in position the complete the 

Proof of Proposition 13.11 Let y be a potential which coincides with V on B{ai, \xq) and 
such that the equivalent of ([9]) for V holds everywhere on M'^. Consider the unique solution 

Ve of 

dtVe-d,,Ve = -\vV{Ve) (3.7) 

on A, which coincides with Vs on dpA. By scalar multiplication of (j3.7p with Vs — en and 
taking into account ([9]), we are led to 

on A. It follows from the maximum principle that \v£ — ai\^ achieves its maximal value on 
9pA, and hence that E -B(o"j, on the whole A. Since V = V ois. B{ai, \xq), we deduce 
that Vg = Vs on A, which is the desired conclusion. □ 
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4 Estimates off the front set 

In the previous section, we have proved Proposition 13.11 which corresponds to the first part 
of the claim (namely ()29p ) in Proposition [H The purpose of this section is to present a result 
which will yield the second part of the claim, namely (j30p . 



Proposition 4.1. Let Vg be a solution to (jl7p verifying assumption {Hq). Let x E M, t > 0, 
r > and s > t such that 

Veiy,T) e B{ai,\io) for all {y,T) e [x - 3r/4,x + 3r/4] x [t,s]. 

Then 



[ ^ ' ee(i;e(y,s))d2/ <Mo(l + 4^ 

Jx-r/2 \ 



exp(-^(.-t) + 2i^-^exp(-W^--) 



Proof. Let < 93 < 1 be a smooth cut-off function such thaJll ip = outside [x— 3r/4, a;+3r/4], 

ip = 1 on [x — 2r/3,x + 2r/3], \ip'\ < Qir/r and \{p"\ < 727r^/r^. Let Wg = dxVs, so that Ws 
satisfies the equation 

dtWe - d^^W, = -\v^V{Ve)w,. (4.1) 

We sellfl fe = ifWe , so that 

dtfe - d^Je + \v^V{Ve)fe = -1^' d^W, - p>" W, (4.2) 

on M X [t, s]. By scalar multiplication with 2/^ we are led to 



dt\fe? - + 4 W(i;,)/e • /. = -1\dxfs? - A^'f'wAw, - 2^"^\w,\\ (4.3) 



Since (fdxW^ = —ip'we + dxfe, we have 

^^\dxWs\^ < 2\dxfe\^ + 2iip'f\We\^ 

and since 

-Aiptp'wedxWe < ip'^\dxWe\'^ + 4{ip')'^\Ws\'^ , 

we obtain 

dt\fe\^ - dxxlfel^ + ^V^V{Ve)fe • /. < {di^'f + 2\^"\p) \We\^ . 

In view of the pointwise bounds on p and its derivatives, as well as of Q, we are led to 

dt\fe\^ - dxx\fe\^ + ^V^V{Ve)fe ■ fe < 3607r^r-^\We\^lA, (4.4) 

where A = [x — 3r/4, x — 2r/3] U [x + 2r/3, x + 3r/4] contains the support of (p' . 



^The constants hereafter were computed for a sine transition function. 
^Notice that is implicitly defined on the whole R x [t, s]. 
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It follows from the comparison principle and Duhamel's formula that 
\f,\\.,s)<expi-^is-t))Gs^t*\MH;t) 



+ SeOvrV-^ ^ exp(-^(s - t))G,_, * \we\hA{;T) dr, 



(4.5) 



where 



In particular, it follows from the triangle inequality that 

A7 



||/e(-,^)lli2(s)<exp(--V(s-t))||/.(-,t)|| 



+ 360^V-2 f exp(-%(s-r))||G,_,*|^^.l'U(-,r)||ii(B)dT, 
Jt £ 



(4.6) 



where B = [x — r/2, x + r/2]. 

Since d{A,B) > r/6, we have, for t < t < s, 

\\Gs~T * \Ws\'^lAi-,T)\\L^B) < \B\\\Gs~T * |We|^lA(-,r)||ioo(B) 

< \B\ma^\\Gs-r\\L°°{A)\\\w£\^'^A{-,T)\\L^B) (47) 

where we have used the convolution inequality ||(7 * /i||oo < llsHooH^Hi and assumption (i^o)- 
Combining (j4.6p with (j4.7p we obtain 



.(•,s)||i2(s) <2e-iMoexp(-^(s-t)) 



+ 360^lr-2e-iMo exp (-^(^ - r) - ^4(7=7)) 7= 



We next use the inequality x < exp(x /288) in order to bound 



[ 144(s~r) J v^537 ^ V 288{s-t) J ' 



and then the inequality 



in order to obtain 



7^^"" " + 288(s - r) - V Tei 



A, 



exp ( --^(5 - r) - ) ^ < exp(-W -f- 



X7 r 
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Inserting the last inequality in (14.81) , multiplying by | and integrating yields therefore 



x+r/2 1^ „, 12 

e^-^{y,s)dy<Mo 

x-r/2 ^ 



S - t 



A- r 



exp(-^(.-t))+2^^— Iexp|-y:^- 



(4.9) 



A completely similar computation for the function = (p{vs — Ci), in view of ([9]) yields the 
estimate 



x+r/2 Y( \ \ + 

^ ''{y,s)dy<A-^Mo 



„ -r/2 e ~" ■ " A, 
The conclusion follows by summation. 



exp {-^{s- t)) + 2^^^^ exp 



Aj ^ 
2 6e 



. (4.10) 
□ 



Proof of Proposition [T] completed. It is an immediate consequence of Proposition 13.11 
and Proposition 14. li □ 



5 The notion of ^-confinement and optimal coverings 

Definition 1. Let X be a metric space and S C X. Given a collection of distinct points J 
in X , p > and < k < 1, we say that {J,p) is a n-confined covering of S at scale p if 



and 



S C UaejB{a,Kp) , S n B{a, Kp) ^ $ M a £ J 



dist(a, h) > K, ^ p Ma in J 



(5.1) 



(5.2) 



We would like to draw the attention of the reader to the fact that Definition [T] involves in 
particular two parameters: the first one, p is the typical length scale at which the confinement 
takes place, whereas the second one, k, controls the rate of the confinement. 

A simple combinatorial argument yields 

Lemma 5.1. Let X he a metric space, and consider i distinct points ai,...,ai in X. Let 

6 > and < k < 1 be given. Then there exists p > such that 



6<p< 



(5.3) 



and a subset J of {ai}i<i<i such that {J, p) is a K-confined covering of S = {oi,--- , a^} at 
scale p. 

Proof. The proof is by iteration of merges in a finite number of steps. First, consider the 
collection J = {ai, • • • Obviously, (jS.ip is satisfied whatever our choice of p > is. If 

(15. 2p is verified with the choice p = 5 then there is nothing else to do. Otherwise, there are 
two points, say oi and 02 after a possible relabelling, such that 



dist(ai, 02) < K. ^5. 



(5.4) 



We then consider the collection J = {02, . . . ,ai] and set p = k ^J. Since by assumption 

dsaD 

ai G B{a2, up), 
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condition (15. ip is therefore verified. As above, either ()5.2p is verified for this choice of p and 
J or we go on in the same way. If the process does not stop in £ — 1 steps, at the — 1)"^ 
step we are left with one single element in J and p = k^'^^^"'^'^ 5. At that point ()5.2p becomes 
void and hence satisfied. □ 

We now apply the previous result to the front set. An easy consequence of Corollary [J and 
Lemma 15.11 is 

Lemma 5.2. Let e > and u he such that 

£e{u) < Mo < +00. 

Given 6 > and < k < 1 verifying k5 > 2e, there exist 

S<p< (|)-2M>^o^ 

and a finite subset J C T>{u) such that [J(J) < Alo/rjo and {J,p) is a K-confined covering at 
scale p ofV{u). 

Proof. By Corollary [T] and (fTSj) we deduce that there exist I < ^ and £ points xi,...,Xi 
such that T>{u) C yjl^i[xi — e,Xj + e]. Lemma [5TT] used with ^ yields a ^-confined covering 

(J,p) of {xi,...,X(] at scale p, with 5 < p < "o since 'iJ < I < Therefore, 

I'(ii) C \Jxj<^j[xj — h,Xj+h], where /i = f/O+e < k/O. Since moreover > 2k^^p > K^^p, 

{J, p) yields a K-confined covering at scale p for V{u). □ 

Definition 2. For e, tt, (5 and n as in LemmalKE, we set 

n{u,K,5) = inf{tj(J)}, 

where the infimum ranges over the sets J for which {J, p) is a K-confined covering of T){u) 
for some 5 < p< (f )-2*^o/^o5. 

6 Stopping times 

In this section we define a notion of exit time for the front set, state and establish a result 
which is the core of the proof to Theorem [TJ 

In the whole section, denotes a solution to ()17p verifying (Hq). We fix a time t > and 
a length scale 5 > 0. We also fix the value of k to 

o -1 . 1 

° 16Afo ~ 16 

We assume that hqS > 2e, so that Lemma [5]2] yields the existence of p such that 

and J = {ai,--- , a^} C T>{t)) such that (J, p) is a KQ-confined covering of V{t)) at scale 
p. Without loss of generality, we may assume that this covering is optimal in the sense of 
Definition^ i.e. j:)(J) = n{v£{-,t), k,q,S). 
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6.1 Defining stopping and exit times 

We define the exit time Ti G [i, +00] by 

Ti = Ti{vs,t,{J,p)) =mi{s>t s.t. V{s) ^UaejB{a,p)}, 
and the dissipation time T2 E [t, +00] by 

T2 = T2{vs,t) = inf{s > t s.t. e / \drVs\'^{x,T) dxclr > 

Jt Jr 8 

Finally, we define the target time T3 by 



2 1 ______ kvp, 



where 



r3.r3(t,e,,,M.)=t + /^^=_exp(^J), 



= 215(1 + 4max^), = mm{min(A7 , W ^1)} , 

4=1 A- «=i 1/2 6 



6.2 Dissipation or splitting 

QQ 

Set Po = 4ao and 

70 = max{ao/3o, -j^ ^og{4:alKv), ^/GKyao/kv}- 

The main result of this section is 
Proposition 6.1. Assume that 6 > 706. IfTi < T2 and Ti < T^, then 

n(ve(:,T{),Ko, — ) > n{ve{-,t),Ko,6) + 1, (6.1) 
Po 

where = Ti — AnQaQ^^d'^ > t. 

In other words, if dissipation does not occur for a sufficiently long time, a front may exit 
the confinement if and only if has already split before to yield an additional (well separated 
at a finer scale) front. The small shift in time in the definition of is motivated only by 
technical reasons. 

Proof of Proposition 16.11 We divide the proof in a number of steps. 
Step 1: lower bound on the exit time. Set = t + AKQaQ^p"^. We have 

Tj^ > , or equivalently, Ti > t + SkqUq ^p'^ . 

Proof. By definition of Ti, it suffices to prove that for t < s < t + SkqQq'^p^, 'E'(s) C 
Ua£jB{a, p). If X ^ UaejB{a, p), then since 'D{t) C UaejB{a, kqp) and kq < jq, we have 
[x — 4:K,qp, X + 4:Kqp] n T>{t) = 0. It follows from Proposition [T] that x ^ 'P(s) for t < s < 
t + IGkIoq^ p'^ . □ 
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Step 2: localising an exit point. There exist 6* G "Di^Ti) and G {1, • • • ,i} such that 

dist(6,,aij G + 3ko)p]. (6.2) 

Proof. Let x G R \ [JaejB{a, p + Skqp). Since T^" < Ti, it follows that [x — 3kq6, x + SkqS] n 
T>{T{) = 0. We induce from Proposition [U that x ^ T>{s) for t < s < Ti + bK^aQ^p"^. By 
definition of Ti, this imphes that there exist a sequence of times (ij)jgN in [Ti,Ti + 5KQaQ^5'^] 
such that tj decreases to Ti as j ^ +oo and a sequence of points (6j)jgM such that bj G 
and dist(6j, U^^iloj}) G [p, (1 + 3kq)p] for each j G N. Let 6* be an accumulation point of the 
sequence {bj)j^n. Then G V{Ti) and there exists i* G {1,...,^} such that dist(oi^ , 6=,,) = 
dist(6, Uf^i {a-}) G [p, (1 + 3ko)p]. □ 



Step 2bis: There exist c=k G T^{T^ ) such that 

dist(c*, 6*) < 2kqp. 

In particular, 

2 4 

dist(c=^,aij G [-p, -p\. 



(6.3) 



(6.4) 



Proof. Inequality (|6.3p follows directly from Proposition [T] applied at time Tj^ with r = 2ko5. 
Since kq < 1/16, ([631) then follows from {62]) and ([O]). □ 

Step 3: confinement of the energy after a boundary layer in time. Let x G M \ 

UagjS(a, 3kop)) then 



X + KQp 
X — KQP 



esivsiy,t^))dy < KvMoex.p{-kva(^ -), 



where t+ = t + 4kqQq p- 

Proof. Since P(t) C [JaejB{a, nop), it follows that [x — 2ko/5, x + 2ko/o] H P(t) 
Proposition [1] with r = 2ko/0 we therefore obtain at s = 



X + KOp \f 

ee{ve{y,t+))dy<Moil + A^ 

X-KOp -^j 



exp(-ATao-« ^ ' 



+ 2^^ao^exp(- 



(6.5) 
• Applying 

>^J_Kop^ 

2 3e ^ 




where j G {1, • • • ,(?} is such that Ve{x,t) G i?(cTj, (J-o)- The conclusion follows from the 
definition of Ky and ky noticing that since kq5 > 2aQe, we have 

— ^ > 4ao . 



□ 

Step 4: energy estimates off the front set. For t < s < Ti and for z G {1, • • • , ^} we 

have 



ee{ve{y, s)) dy < KyMo exp{-ky — 5-) H 5- exp(-A;v-^ 



(6.6) 



The same estimate holds for integration on the interval [oj — 2p, Oj — ip]. 
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In case < , we have therefore c^, < ~\p (!6-4p , and this imphes 



Proof. It suffices to apply Proposition 14.11 with x = ai + 2p and r = |p. Indeed, we have 
[x — |r, X + |r] = [ai + /?, + 3p], and this set is disjoint from V[s) whenever t < s <Ti,hy 
definition of Ti. □ 

The next step contains the core of the argument. 

Step 5: existence of a splitting. There exists c* G 'D{T^) n B{ai^,3p) such that 

dist(c„c*) > |. (6.7) 

Proof. Assume by contradiction that there are no such c*. We wiU obtain a contradiction 
using the evolution equation (j24p for a suitable test function x- We distinguish two cases. If 
c* > fli, , we have c* > a-i^ + |p by ()6.4|) . and this implies that 

P(Tf ) n [ai, - 2p - 2koP, fli. + 2ko/9] = 0. 

We may thus apply Proposition [T] exactly as in the proof of Step 3, obtaining 

/ e,(7;,(y,ri))dy < E:yMoexp(-%ao -) (6-8) 
for every x S [aj, — 2/?, aj,]. 

therefore c,,, < — 3 
P(Tf ) n [a,, - 2koP, fli. + 2p + 2kop] = 0, 
so that we obtain as before 

/ e^iveiy, Ti)) dy < KvMo exp{-kva^^l^) (6.9) 

Jx — KQP ^ 

for every x G [oj, , Oj, + 2/?] . 

Since 6* G 'D{Ti), we have by the clearing-out lemma, 

e,(7;,(2/,Ti))dy >7?o. (6.10) 

Let X be a smooth function with compact support in [ai^ — 2p, ai, + 2p\ such that xiv) = 
y-au for ?/ G [ai,-|p,ai.+|p], x(?/)(y-ajj > for y G M, ancil||x||oo < 2p, \\x"\\oo < 24p-i. 

We integrate equality ()24p between t"*" and Ti and estimate each of its terms. 

Combining (|6.8p or (|6.9p with (j6.10p . using the fact that by (|6.2p \ai, — y\ > ^p for y £ 
[6* — £,K + e], and in view of the properties of x, we are led to 

1/ es{ve{y,Ti))x{y)dy\ >l-r]op-2KQ^KvMoexp{-kvaQ^-)p. (6.11) 

We now estimate the same integral at time t~^. Combining (j6.5p outside L)aejB{a,3KQp) 
with the bound (-ffg) inside Ua£jB{a, 3kqp), we are led to 



/ 

Jm 



,{v,{y, t+))x{y) dy\ < Zn^M^p + 2k^^KvM^ exp(-%ao 3^)p. (6.12) 



These last estimates are fulfilled for a (regularization of a) degree three interpolation polynomial. 
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r 



We also have, since ||x||oo < 2/?, 

I I e\dsVe{y,s)W{y)dyds\ < < -riop, (6.13) 

where we have used the fact that by assumption Ti < T2. 

Finally, since ||x"l|oo < 24/5"-^, since x" is supported in [aj^ —2/3, — |/3]U [cj, + Oj, +2/?], 
and since |'?£(a;)| < pointwise, we obtain from (16.61) . 

Jr V e p e / 

provided t < s < Ti. Integrating the last inequality from to Ti we are led to 

Uveiy, s))x"{y) dyds\ < 48p~'KvMo + eM-kv-)) • (6.14) 

Combining (f6TT]) . ([6A2]) . ([6l^ and ([6l^ with ([Ml) we deduce 

- 3KoAfo - - 4^0^^!/ Mo exp(-A;yag ^^)p 

< \%p-^KyM^ (f- + exp(-fcy^)) . (6.15) 

\kv 2p^ £ J 

Since by assumption ^ > ^y GKyao/kv , we have 

3 

Since by assumption | > log(4aQi^'y), we also have 

4kq^ KyMo eM-kvao''^)p < '-^p. (6.17) 

e 64 

Combining ()6.15p with ()6.16|) and ()6.17p we finally deduce 

A8p-^KvMo^^^^ffieM-kv^) > ^P, 

£ 62 



so that 

Ti>t + p \ J- — ^M^--) > T3, 

2y/bKvao 2 £ 

the desired conclusion. □ 
Step 6: persistence of the other fronts. For i G {1, • • • \ 

V{T^)n[ai-2p,ai + 2p]^^. 
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Proof. The proof is by contradiction, very similar and actually simpler than the one of Step 
5, so that we only briefly sketch it. It relies also on equation ()24p . with a function x which 
is here taken to be non negative, identically equal to 1 in [oj — |p, + with compact 
support in [oj — 2/3, Oj + 2p] and verifying ||x||oo ^ 1 and ||x"||oo 

< 145-2. This time, the 

dominant term is given by / e^x taken at the initial time t, and not at the exit time. □ 

Step 7: proof of (|6.ip . In view of Step 5 and Step 6, there exists i + 1 points in DiT^) 
such that the mutual distance between any two of them is bounded from below by |. Any 
covering of T>(T^) with balls of radius smaller than p/A must therefore contain at least i + 1 
balls. Notice that the existence of a KQ-confined covering of T>{T^) at scale p' < p/A follows 
by Lemma 15.21 with 5 replaced by 5' = S/Pq: by assumption we have ko6' > 2e and the scale 
p' of the KQ-confined covering does not exceed J'(M-2*^o/''o) < 5/4 < p/4. The conclusion 
follows. □ 



7 Proof of Theorem Ibis 

We distinguish between two cases, depending on the value of the ratio 

Case I : i? is sufficiently largj^ with respect to e. 

This is the most interesting case, where we apply the analysis of the previous sections. In 
Section [6l the initial time t as well as the typical length scale 5 were fixed. In order to prove 
Theorem Ibis in the present case, we will use Proposition 16.11 inside an iteration argument, 
containing a finite and bounded from above number of steps, each step corresponding to a 
new time tn and a new length scale 

More precisely, we set to = and will fix the value of 5q at the end of the argument. It 
remains to define the times t„, and the scales 5n iteratively. 

Let n > be fixed. In view of Section [H provided 6n > 7o£ there exist 

Sn < Pn < (|^)^^'^»/*5n (7.1) 

and Jn = {a",-- - , } such that {Jn,Pn) is an optimal KQ-covering of V{tn) at scale pn- 
Define the corresponding stopping and exit times Ti^^, T2^rn ?3,n according to Section [6Tl 

We distinguish three cases : 

(Case 1)^ : T^^n < 2i,ti In that case, we define tn+i = ifin = ^3,n and the iteration 
process stops. 

(Case 2)^ : T^^n > 2i,ti and T2,n < T\^n In that case, corresponding to a dissipation 
time, we set = T2^n and 5n+i = (we reset the resolution scale to its initial value). 

(Case 3)^ : Ti^„ < min(T2,Ti5 ^3,^) In that case, corresponding to a splitting, we set 

tn+l = T^^n 5n+l = J^^n- 

By construction, in each case, 

V{s) C V{tn) + [-2p„, 2pn] for tn<S< tn+l. (7.2) 
^''The actual limiting ratio is given by (|7.5p . 
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Moreover, if (Case 3),„ holds, then by Proposition 16.11 

n{Ve{-,tn+l),K0,6n+l) > n^Vei' , tn) , Kq , 6n) + 1- (7.3) 

Since 6n+i > e by assumption, n{v£{-,tn+i), Ko,Sn+i) < ^ Corollary m It follows 



by (|7.3p that (Case 3)„ may only occur for at most consecutive values of n, provided 5„ 
remains greater than jqe. In view of the iteration process, this also implies that 5n is bounded 
from below by 

Sn > Po '^0. (7.4) 

On the other hand, since the total dissipation is bounded from above by the total energy 
of the initial datum, (Case 2)„, in view of the definition of T2^n, may only occur for at most 

4— distinct values of n. 

vo 

We induce from the two previous limitations that the process may contain at most 

+ l)(4^*) + 1 

Vo Vo 

steps, so that t^^ is weh defined. Moreover by ()7.2p . (|7.ip . the fact that ^'kLof^o'' — 2, and 
the triangular inequality, we obtain 

V{s) C V{0) + [-rfin, rfin] for < s < tf,n, 

where 

2 \ Vo J 2 



We may now define the value of Sq to be 

an 1 
6o = Ra, «(^ + 4)-i 

If 



i?>(ao)«(Y + 4)/?o''°7oe, (7.5) 

we infer from (|7.4p that the condition 6n > 7o£ is met along the whole process. We also infer 
from (j7.4p and the definition of T-^^n that 



where where 



tfi„>(^j exp(--) 



Ko = /^o"" af (^ + 4) max{ ^24Kyao, 7^} 

2 Kv 

/ ,^ N 2{^)(1+^) ^ 

> 2^° 32— max{max( 
V m J *=i 

The proof is then completed in case (I7.5|) holds. 

Case II : (17. Sp does not hold but R > aos. 

In that case, it suffices to invoke Proposition [TJ The conclusion follows adjusting the 
constant Kq if necessary. □ 
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8 Proof of Theorem [2] 

It is a direct consequence of Theorem Ibis and Proposition l4.1l that for < t < (■^)^ ^^p(t^^ 
and X ^ [xq — ^R, xq + ^R], 



x+e 

es{v,{y,t)dy<KMo 

x—e 



for some constant K depending only on V. To derive the pointwise bounds given by ()2ip . it 
suffices then to invoke parabolic regularization and scahng. Indeed, if u is a sohition of ([1]) on 
a cyhnder of the type [y — 1, y + 1] x [s — 1, s], such that v € B{ai, \xq) for some i £ {!,• • • ,q} 
on that cyhnder, then 

\dtv{y,s)\^ + e{v{y,s)) <K' [ e{v{z,T)) dzdr. 

Therefore, by scahng if -Wg is a solution of (jl7p on a cylinder of the type [y — e, y + e] x [s — e^, s], 
such that Ve G -B(o"i, |J.o) for some i G {1, • • • , on that cylinder, then 

£^\dtVe{y,s)\''+£ee{Ve{y,s))<K'£~^ j r^\,{Veiz,r))dzdT. (8.2) 

J s—e^ J y—e 

The conclusion with Ki = 2Km.ayi{K\ 1) follows combining ([8T]) and ([812]) . □ 

9 Relaxation to stationary fronts 

The aim of this section is to provide a proof to Theorem [3l The starting idea is to determine 
a good time slice for which the integral of the dissipation |9jUep is small. Then, the main 
part of the proof is devoted to the study of solutions to the perturbed ordinary differential 
equation 

Uxx = e^'^VV{u) + f onM, (9.1) 
where the function / belongs to L^(]R). 

9.1 Study of the perturbed equation ( 19.11) : initial value problem 



It is useful to recast equation (19. ip as a system of two differential equations of first order. 
For that purpose, we set w = eux so that ()9.ip is equivalent to the system 



1 

Ux = -W 

I (9.2) 

Wx = -VV{u) + £f, 



which we may write in the condensed form 



Ux = ^G{U) + eF on M, (9.3) 
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where, for x in M, we have set U{x) = {u{x),w{x)) and F{x) = (0, /(x)), and where G 
denotes the vector field on M.'^^ given by G(«i,U2) = {u2,'VV{ui)). Notice that 

\VG{ui,U2)\ <^(|^xi|), 

where A > 1 is some continuous non-decreasing scalar function. We next compare a given 
global bounded solution u of ()9.ip to a solution of the unperturbed equation 

u,, = e~^VViu) , (9.4) 

with similar initial condition at some point xq £ M. We denote accordingly = {u^,e~^u^) 
on its maximal interval of existence. 

As a consequence of Gronwall's identity, we have 
Lemma 9.1. Let u and uq be as above. Assume that for some xq £M and a > 0, 
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- U\0)\ + f Wfh < exp f_ Mn\\oo + l)a \ ^^^^^ 

V2^(||u||oo + 1) V e / 

Then is well defined on [xq — a, xq + a] and we have 

l|f/-f/»ll.^„«-«.,, < (l^(0) - U'ml + exp (4*^ 

(9.6) 

Proof. Without loss of generality, we may assume that xq = 0. Let / be the largest interval 
containing and such that 

< ||ii||oo + 1- (9-7) 
On /, since {U - U^)^ = G{U) - G{U^) + eF we obtain the inequality 

It follows from Gronwall's inequality, that, for x £ I, 
\{U - C/0)(x)| < exp ( -4(ll^lloo + l)kh _ f^O)(o^| 



- y)| exp ( \ dy\, 



+ 

so that by the Cauchy-Schwarz inequality, we are led to the bound, for x £ I, 
\m - U')ix)\\ < (mO) - C/°(0)| + f Wfh] exp Mdlnll +l)|x| \ 

Hence, if ()9.5p is verified, then [—a, a] C I and ()9.6p follows. □ 
We will combine the previous lemma with 
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Lemma 9.2. Let u be a global solution of ()9.ip such that £eiu) < Mq < +00. Then 

Ueiu)\\oo < V2eMo\\fh. 
Proof. This is a direct consequence of the equaUty 

^Un) = ef-u, 

Cauchy-Schwarz inequahty, and the fact that £,e{u) tends to zero at infinity since u has finite 
energy. □ 

Corollary 9.1. Let u be a global solution of (j9.ip such that £e{u) < Mq < +00. There exist 
a constant Co > depending only on Mq and V such that if xq E T^iu) and if 

b = u TT log f T I > (9-9) 



then there exists a solution of (j9.4p defined on [xq — b,XQ + b] and verifying 

Uyl") ^ 0, (9.10) 

||t/-C/°|koc([,„_,,,„+,]) <exp(-^i^^^^^(6-G)) , (9.11) 

for every < a < b. 

Proof. Since xq S T>{u), there exist a constant cq > (depending only on the choice of |j.o 
and the eigenvalues ) such that 

e|n,(xo)|V2 = ^^^^^ + Uu){xo) > ^ +e.(^)(^o) 



so that by Lemma [921 

e|n.(xo)|V2> j-V2eMo||/||2. 



Since a > 0, we infer that CqE^ \\f\\2 < 1 and therefore 

e|n,(xo)|V2 > ^ - V2eMoC^^e-h 

We first require Co > 2^/2Mo/co so that we obtain, since < e < 1, 

£\u^{xo)\ > V^. (9.12) 

We let n'^(xo) = u{xq) and we wish to define u^{xo) in such a way that S,eiU'^){xo) = 0. This 
may be achieved in general in a non unique way. We choose u^{xo) as the unique positive 
multiple of Ux{xq). It follows from the equality 

e\u^ixo)\'^ - e\ul{xo)\'^ = 2Cs{u)ixo), 
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from the bound ()9.12p . and from Lemma 19.21 that 



\e{u,-ul) {xo)\<2V2^e^\\f\\2, 



or, since u{xq) = u^{xo), that 



\U{xo) - U\xo)\ < 2V2^e^f\\2. (9.13) 



In order to apply Lemma |9.H we estimate 
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1^(0) - U'{0)\ + ^=^^—=11/112 < ( 2^2. ^-^^ - 



V2^(||n|U + l)"" V"'"V^ ' 2 
We next require that Co > 2\/2^ + ^ so that by definition of b, 

The conclusion then follows from Lemma |9.1[ □ 

9.2 Study of the perturbed equation (19.11) : boundary value problem 

In this subsection, which may be viewed as the elliptic counterpart of Section HI we prove 
estimates regarding a solution u of (19.11) on some interval, provided that this interval has an 
empty intersection with the front set of u. 
More precisely, we have 

Lemma 9.3. Let u be a solution to (|9.ip such that £e{u) < Mq < +oo, and let xq £ M, r > 
be such that 

[xq — r, xq + r] n V^u) = 0. 
There exist a constant Ci > depending only on Mq and V such that 

eesiu) < Ci (ei||/||2 + ^exp(-^)) (9.14) 

on [xq — r/2, xq + r/2]. 

Proof. The proof is very similar to the one of Proposition 14.1^ so that we skip part of the 
details. In the sequel, C denotes a constant which depends only on V and whose actual value 
may vary from line to line. If i E {1, • • • , g} is such that u{x) G i?(CTi, jJ-o) on [xq — r,xo + r], 
then the function w defined by 

w{x) = {u- aifx^. 

where < x ^ 1 is a smooth cut-off function equal to one on [xq — |r, xq + \r] and with 
compact support in (xq — r, + r), satisfy on M the differential inequality 



^w^Lvj^or isupp(x')- 



-w^xx + <Cf + Cr^H 
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By the comparison principle, we obtain 



w < C{K *f + r '^K* lsupp{x' 



where 



We then estimate 



e 

K{x) = — ■j= exp 



\\K*f\\^<\\KMfh<Ce-^\\fU 
and since dist([3;o — r/2, xq + r/2], supp(x')) ^ 

Ik"^^* lsupp(x')lli°°([^-o-§,^o+5]) ^ ^"^ll^llL-(R\[-f,5])meas(supp(x )) 



< C- exp — V 

r V * '^^ e 



Since x = 1 on [xq — r/2, + r/2], we therefore obtain, for x £ [xq — r/2, xq + r/2], 

V{u{x)) < C (.§ Il/lb + exp(-^)) . (9.15) 
Finally, by definition of we have 



eUu) + V{u) 



2 

so that for x G [xq — r/2, xq + r/2], by (19.15P and Lemma W?2\ 

e'\u,\\x) < V2Moe2||/||2 + C (et H/lb + exp(-^)) ■ (9.16) 

Combining (j9.15|) . (|9.16|) and the fact that e < 1, the conclusion follows for the choice 
Ci = C + y/2Mo. □ 

9.3 Proof of Theorem d 

Since the total dissipation of energy is bounded from above by Mq, by averaging, there exist 

< T < (— )2exp 



such that 

e jjdtVe{x,T)\^dx < KlMoR^^exp {~^) 
Hence, Vir{-,T) satisfies equation (j9.ip where 

b < Ko^JW^R-h-'^ exp (-^^) • 
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In order to apply Corollary l9.lt and in view of ()9.9p . we first estimate 

e . f I \ 

log 



A{\\v,{.,T)\\^ + l) '^"{Coel 



2 , 



- , . I R / R 

— T7T\ — 7 — T ' — 7= — ^■'^P 



Ai\\ve{-,T)\\^ + l) "VCoi^'oV^e \2Koe 



(9.17) 



>2A 



provided K2 is large enough depending only on Mq and V. Next, we apply Lemma 15.21 with 
K = J and 6 = ^2 to . This yields 

R _4Ho R 

— 2 TO < r < 

K2 ~ ~ K2 

and a collection of points {aj}j^j such that {jj < AIq/tjo, 

V{T)cUjej{aj-r/A,aj+r/A) 

and 

dist(ai, flj) > 4r \/ i ^ j £ J. 

In particular, 2r <b and for each j £ J we may apply Corollarv 19.11 with j;o = and a = r. 
This yields, after rescaling, the collection {Uj}j^j and the estimates 

||T;,(-,r) - U,{'-^)\\+e\\d. (^Ve{;T) - U/-^)^ \\ < exp [-A{\\vs\\oo + 1)J) 

< K2 exp 
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in the space L°°{[aj — r, aj + r]), provided K2 is large enough, depending only on Mq and V. 

Consider now xq G M such that dist(x, Ujg j{aj}) > r. Since [xq — jr,xo + |r] nP(r) = 
by construction, we obtain by Lemma 19.31 



eesive{xo,T)) < Ci ( Kq^/Mq- exp -7-^ + ^ exp 



so that 

||7;e(2;o,r) - OjW + £\\dx{vs{xQ,T)\\ < exp 



provided K2 is large enough, depending only on Mq and V. The conclusion follows. □ 
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